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Starting with a time-independent Hamiltonian h and an appropriately chosen solution of the 
von Neumann equation ip(t) = [h, p(t)[ we construct its binary-Darboux partner h\(t) and an exact 
scattering solution of ipi(t) = [hi (t), pi(t)\ where hi(t) is time-dependent and not isospectral to h. 
The method is analogous to supersymmetric quantum mechanics but is based on a different version 
of a Darboux transformation. We illustrate the technique by the example where h corresponds to 
a 1-D harmonic oscillator. The resulting hi(t) represents a scattering of a soliton-like pulse on a 
three-level system. 
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One of the main ideas of supersymmetric (SUSY) 
quantum mechanics (QM) can be summarized as follows 
0. Assume we know a ground state 1-00) of a stationary 
Schrodinger equation (SE) 

H\iPo) = (H km + V)\ifo) = E \ip Q ) (1) 

with some V and Eq. Using |^q) we construct an "anni- 
hilation" operator A = A(ipo) satisfying A\ipo) = and 
H - E = A^A. Now define |-0i> := A\ip) (here \ip) is 
any eigenstate of H linearly independent of |"0o) j with 
eigenvalue E) and Hi = AA^ = H\^ n + Vi. Hi is the 
so-called SUSY partner Hamiltonian of H. Then, using 
AH = AA^A = H\A, one finds that 

(ffkin + Wl) (2) 

In a single step we have produced a new potential V\ and 
one solution of the corresponding stationary SE. 

The map V — > V\ is known to be particular example 
of a Darboux transformation (DT) . All DT transform 
a "potential" V into V\ and simultaneously generate an 
"annihilation" operator A(ip ) satisfying A(ip )i(;Q = 0, 
where ipo is a solution of some partial differential linear 
equation associated with V. The physical interpretation 
of such an abstract "potential" depends on the problem. 

SUSY QM deals with linear SE and for this reason the 
density matrix generalization is not interesting: Hi can 
be inserted either into the SE or into the von Neumann 
equation (vNE) ip = [Hi , p] . 

However, the vNE has a structure which is alge- 
braically different from this of the SE and therefore allows 
for different DT. A candidate is the so-called binary DT 
(BDT) originally constructed in ||] and applied to opti- 
cal soliton equations. Quite recently the technique was 
applied to Yang-Mills equations Q] and nonlinear vNE 
|^],[|. A tutorial introduction to density matrix applica- 
tions is given in fffl. There are formal analogies between 
the BDT and the "dressing method" of Zakharov et al. 



P| but technically the two procedures are inequivalent 
(for a discussion cf. 

The purpose of the Letter is to show that the BDT 
leads to a new kind of SUSY-type QM for density ma- 
trices which does not have a counterpart in SUSY QM 
based on SE. At an intermediate stage of the construction 
we solve the nonlinear vNE 

ip=[H,p% (3) 

where H is a time-independent Hamiltonian. The set 
of solutions of (j^) contains all the pure states of stan- 
dard QM since for p 2 = p (||) reduces to the linear 
vNE. For p 2 7^ p there exist at least two more classes 
of solutions. One of them occurs for p's satisfying cither 
p 2 — ap = with a € -R, a ^ 1, or a weaker condition 
[H, p 2 — ap] = (now a = 1 is acceptable) . In both cases 
pit) = e~ laHt p(0)e laHt . The second class is of the form 
p(t) = e- iaHt p iat (t)e iaHt with p^t (-00) ? p^t (+00). 
These additional solutions, called here the self-scattering 
(SS) ones, are fundamental to our construction because 
of the following property: Each SS solution of the non- 
linear vNE (H) with the time-independent H is simul- 
taneously a scattering solution of a linear vNE with a 
time- dependent Hamiltonian hi(t). Both the SS solution 
and the new Hamiltonian are algebraically constructed 
in terms of BDT. The construction does not make use 
of supercharges and for this reason the resulting partner 
Hamiltonians will be termed the binary-Darboux (BD) 
partners. 

The BDT method of solving (§) was described in ||U • 
We start with the family of Lax pairs, parametrized by 

lo e C, 

Zu,\il>u>) = (p-uH)\ip ul ), (4) 
i|&»>= (Hp + pH-LjH 2 ^^), (5) 
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where z u is a complex eigenvalue. The pair (H)-(0) is here 
the analogue of (|l]) ; p and h = Hp + pH play the role of 
the "potentials". 

The connection of (f§)-(||) with (||) is two-fold. First, 
the necessary condition for the existence of \ipj) is given 
by (0). Now assume = iVvM) is an y solution of 
(H) - ® with u> = p and some z M . Denote by P M the pro- 
jector on IV^i) and let A S C be another parameter, and 
p any solution of (||). Defining 



Px = 1 



P _ 



P^l + ^-i^) =:C/^C/t (6) 



we find (cf. |,|) 

z A |VA,i) = (pi-A/f)|^A,i) 
# A ,i) = (#>i + Pi# - Ail 2 ) |Va,i> 



(7) 



(8) 
(9) 



The "Hamiltonians" p—\H and p\ — \H possess the same 
eigenvalue za and their eigenvectors are related by the 
"annihilation operator" A (note that ^4(^ M )|-0 M ) = 0). 
However, these are not the physical BD partners we 
are interested in. BDT transforms the two "potentials" 
p — > pi, h — * hi in such a way that 



ipi 



[Hpi +pii?,pi] = [hi, pi]. 



(10) 



since this condition has to be satisfied whenever \ipx,i) ex- 
ists. The BD-transformed Lax pair can be used 
to repeat the procedure: pi — * p2, hi — ► /i2- 

To explicitly show that the construction of hi is non- 
trivial we have to make an assumption about the Hamil- 
tonian H . We shall concentrate on the isospectral family 
of the 1-D harmonic oscillator (HO) since for Hamiltoni- 
ans with equally-spaced spectrum a strategy leading to 
nontrivial solutions was worked out in detail in [^J. An 
alternative strategy was described in j(| and applied to a 
concrete example in 0. In both cases the result is a SS 
solution. 

We take the Hamiltonian H = eN, where e is some 
parameter, 



N = VVr + n)\r + n){r + n\ 



(11) 



n=0 



and r e R (e.g. for 1-D HO r = 1/2; for 3-D isotropic HO 
r = 3/2). In the Hilbert space spanned by {\r + n)}^L 
consider a 3-D subspace spanned by three subsequent ex- 
cited states \k), |fc + l), and \k + 2). It should be stressed 
that the same strategy can be applied to any H with dis- 
crete spectrum provied there exist three eigenvalues of H 
satisfying E k -Ei = E x - E m . 

In order to obtain a SS solution pi(t) one has to start 
with an appropriate p(t). The problem of how to select 



such a p has been discussed in great detail in ||. The 
fact that (JlJ) does indeed solve (||) with H given by ( |Tl| ) 
can be verified by a straightforward calculation. 

We consider a one-parameter family of solutions, 
parametrized by a € R. Physically the parameter turns 
out to control the scattering process. Mathematically it 
parametrizes an initial condition for the solution of the 
Lax pair (|)-(§). We solve (@)-(|) with 



n -i5Ht 



P (oy 5Ht =: W 5 p(0)Wl 



(12) 



p(0) = -(|fc)(fe| + |fc + 2)(fe + 2| 



s/5 



|* + 1)<* + 1| 



(\k + 2){k\ + \k){k 



(13) 



and p — i/e. For later purposes we have introduced 



the unitary operator W a (t) 



■i HI 



( PL 2[) is a solu- 



tion of (|3j) and therefore the necessary condition for the 
existence of \ipi(t)} is satisfied. 5H in ( |l2|) comes from 
[H, p(0) 2 ] = 5[H, p(0)] and the resulting equalities 



ip=[Hp + pH,p] = b[H,p] = [h,p]. 



(14) 



h = 5H can be regarded as the first element of the pair 
of BD partner Hamiltonians we are going to find. The 
initial condition for (f|)-(||) is 



Mo)) = 



</l+a- 



• \k + l) 



Inserting Pj, which projects on \ipi(t)), into (J6J) with 
p = i/e and normalizing the resulting solution to get 
Tr pi = 1 , we finally get the density matrix 



Pi(t) 



2 

E 

u,v—0 



Pl(t) i+u,i+v I k + u)(k + v\ 



(15) 



where the matrix of coefficients in (tL5|) is 



Pi(t). 



15 + V5 



5 m at) 
at) m 5 



with 



at) = 



a 



(2 + 3« - y/5i) V3 + V5< 
9e 2wo */ 5 + (1 + iy/Ei) a 2 

C(t) ~ 3( e 2^/5 +a 2) ' 



iujQt 



and oj = 10e/(15 + y/E). Writing @ as 



Pi(t) 



PuAt)e^ Nt 



(16) 



(17) 



(18) 



(19) 



one finds, for < \a\ < oo, 



2 



Pint(-00) 



Pint (+Oo) 



15+VE 



15+V5 



5 
5 + V5 

4^/5i q 
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I _ 1 1 4V5i 
\ 3 3 

5 0-3 
5 + V5 



-3 







This shows that ( |l5|) is a SS solution. 

The BD partner h\ occuring in ( jlpj ) is non-unique and 
defined up to an operator commuting with pi . This free- 
dom is useful. Set 



h 1 = {H + ecil)pi + Pl (H + ecil) + ec 2 l 



(20) 



with constant ci, c-i. Denoting cr,- £ = \j){k\, and using 
the above explicit solution we find 

fti(t) = H + Hi(t) (21) 

where H — ujq(N + cil) + ec2l and 

fli(t) = ^%(k + ci + l)o- fe+ i, fe+ i 
vo 

+ y(fc + C! + |)(e(«)^fc,fe+l +f(*K+l,fe) 
+ + Ci + l)(C(t)er fcjfc+ 2 + C(t)crfe+2,fc) 

+ y(fc + C! + §)(e(t)o*+i,*+2 + £(t)cr fe+2)fe+1 ). 

£(t), are essentially the Rabi frequencies. The non- 
uniqueness of 7ij was used again to extend the non- 
perturbed part of (|l]) beyond the 3-D subspace. Our 
construction guarantees that ( |l5|) is a scattering solu- 
tion of the corresponding time-dependent linear vNE 
ipi = [h\(t), pi]. Let us note here that the dynamics of pi 
is related to p(0) by the unitary transformation UiWs. In 
general, taking arbitrary and W a we can alternatively 
define the scattering Hamiltonian as 

hi = iUJJl + aU^HUl (22) 

hi is a nontrivial scattering Hamiltonian provided pi(t) 
is a SS solution of (|). 

( ^jj ) represents a complicated time-dependent three- 
level perturbation of a HO. In order to better under- 
stand the kind of interaction we have produced set a = 1 , 
ci + k + 1 = 0, ec2 = — ojqCi, and d = y3 + \/5(2 + 3i — 
v / 5«)/(2v / 3)- The Hamiltonian now reads 

Lu de luJot 



hi(t) = ui a N - 



10 cosh(w i/5) 
10 cosh(woi/5) 



(cfc,fc+l — CTfe + i ! fc+2) 

(c/c+i.fc — Cfc + 2,fc+i). (23) 



One can think of hi as describing a 1-D HO located at 
x = and interacting with the well-known McCall-Hahn 
"sech" optical soliton jlfj|. Let us recall, however, that 



the result is more general and valid for any H with dis- 
crete spectrum provided the 3-D subspace corresponds 
to three equally spaced eigenvalues. Taking different pa- 
rameters in hi we obtain additional terms reminiscent 
of the "sech-tanh" pulse occuring in inhomogeneously 
broadened three-level media JO]. It is interesting that 
for ci + k + 1 7^ the perturbation Hi(t) contains a 
time-independent term proportional to |/c + l)(fc + l|. Re- 
defining the non-perturbed part by 



H' = oj N + ^=(fc + ci + + l)(k + 1| 
v5 



(24) 



we break the equal spacing of the non-perturbed Hamilto- 
nian, simultaneously detuning the highest and the lowest 
levels from u>o and generating a transition with doubled 
frequency 2ojq. 

A general property of (||) is the fact that (H) n = 
Tr Hp n are integrals of motion for any natural n and 
any solution p lol . In particular, this implies that 
the sum of the perturbed eigenvalues of hi is time- 
independent. The same holds for the average energy 
(E) = Tr hi(t)pi(t). However, the eigenvalues them- 
selves may be time-dependent. For c\ + k + 1 = 0, c 2 — 0, 
the eigenvalues of the restriction of hi to the 3-D sub- 
space are and 




This implies that the BD partners h — 5H and hi are 
not isospectral, a situation that may occur in higher- 
dimensional SUSY. 

The figures illustrate properties of the scattering solu- 
tions. Fig. 1 shows the average position of the 1-D HO 
(x) = ^=Trpi(a + at) as a function of time and a. In 
the asymptotic regions the average is 0. For times where 
(x) « the dynamics is effectively given by 

Pout(t) = e-^ Nt Pint (+oo)e^ m . 

As \a\ grows the moment of SS is shifted towards the 
future. For a — or |a| = 00 there is no scattering since 
Pint becomes time independent. 

The asymptotic probability densities in position space 
g(x,t) = (x\pi(t)\x) are symmetric (implying (x) = 0), 
Fig. 2. Such time-dependent probability distributions 
represent a new type of nonlinear effect. We propose 
to term them the Harzians P2[. 

The above effects can be extended to higher- 
dimensional subspaces. One of the possibilities is related 
to the "weak superposition" principle: For any family of 
solutions {pk} of (J3J) satisfying pkpi — for k ^ I, the 
combination p{t) = ^ fe PkPk(Pkt) is also a solution of (^). 
One can generalize the procedure to many noninteracting 
HO and consideration of systems with degeneracy, such 
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as HO with spin, leads to a nontrivial second iteration 
of BDT: p — > pi — > and ft. — > /ii — ► Another 
possibility is related to the Yang-Mills (YM) case. The 
result of Q shows that a class of YM equations can be 
integrated by BDT. The anti-self-dual YM case is alge- 
braically related to Eulcr- Arnold equations |H| which are 
a particular case of (||) as discussed in || . 

Exactly solvable equations with time dependent 
Hamiltonians are a rarity in quantum mechanics. The 
technique we have described leads to a broad class of 
such equations. The example we have discussed, in spite 
of its simplicity, shows the richness and efficiency of the 
method. The resulting three-level dynamics is highly 
nontrivial and physically interesting. We expect the 
method to prove useful in many branches of quantum 
physics. 
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FIG. 1. (x) as a function of time and the parameter a, 
5 < a < 100, which controls the initial condition. The mo- 
ment of SS moves towards the future (past) as |q| grows (de- 
creases). For |q| = oo (a = 0) SS is shifted to +( — )oo (no 
scattering) . 
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FIG. 2. The Harzian. Probability density in position space 
as a function of time for k — 2.5, a — 5, < t < 20. The 
asymmetry of the probability density is responsible for the 
oscillation of (a;) seen at Fig. 1. 
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FIG. 3. Contour plot of the Harzian from Fig. 2 for 
— 25 < t < 60. The continuous transition between the two 
asymptotic states (with symmetric probability distributions) 
is clearly visible. 
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